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Abstract 

When a quantum system is monitored in continuous time, the result 
of the measurement is a stochastic process. When the output process is 
stationary, at least in the long run, the spectrum of the process can be 
introduced and its properties studied. A typical continual measurement 
for quantum optical systems is the so called homodyne detection. In 
this paper we show how the Heisenberg uncertainty relations give rise to 
characteristic bounds on the possible homodyne spectra and we discuss 
how this is related to the typical quantum phenomenon of squeezing. 

1 Quantum continual measurements 

A big achievement in the 70's-80's was to show that, inside the axiomatic formu- 
lation of quantum mechanics, based on positive operator valued measures and 
instruments, 1,2 a consistent formulation of the theory of measurements contin- 
uous in time (quantum continual measurements) was possible. 2-8 The main 
applications of quantum continual measurements are in the photon detection 
theory in quantum optics (direct, heterodyne, homodyne detection) . 9 ~ 15 A very 
flexible and powerful formulation of continual measurement theory was based 
on stochastic differential equations, of classical type (commuting noises, Ito 
calculus) and of quantum type (non commuting noises, Hudson-Parthasarathy 
equation). 5-16 

In this paper we start by giving a short presentation of continual measure- 
ment theory based on quantum SDE's. We consider only the type of observables 
relevant for the description of homodyne detection and we make the mathemat- 
ical simplification of introducing only bounded operators on the Hilbert space 
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of the quantum system and a finite number of noises. Then, we introduce the 
spectrum of the classical stochastic process which represents the output and we 
study the general properties of the spectra of such classical processes by proving 
characteristic bounds due to the Hcisenbcrg uncertainty principle. Finally, we 
present the case of a two-level atom, where the spectral analysis of the output 
can reveal the phenomenon of squeezing of the fluorescence light, a phenomenon 
related to the Heisenberg uncertainty relations. 

1.1 Hudson Parthasarathy equation 

Let 7i be the system space, the complex separable Hilbert space associated to 
the observed quantum system, which we call system S. Quantum stochastic 
calculus and the Hudson-Parthasarathy equation 17 allow to represent the con- 
tinual measurement process as an interaction of system S with some quantum 
fields combined with an observation in continuous time of these fields. Let us 
start by introducing such fields. We denote by T the Hilbert space associated 
with d boson fields, that is the symmetric Fock space over the "one-particle 
space" L 2 (R + ) <g> C d ~ L 2 (R+;C d ), and we denote by e(/), / £ L 2 (R + ;C d ), 
the coherent vectors, whose components in the 0,1, . . . ,n, . . . particle spaces 
are e(f) := exp (-± ||/|| 2 ) (f , /, (2!)- 1 / 2 / ® /, . . . , (n!)" 1 / 2 /®", . . .). 

Let {z k , k > 1 } be the canonical basis in <C d and for any / £ L 2 (R + ; C d ) let 
us set /fc(i) := (zk\f(t)). We denote by A k (t), A\(t), A k i(t) the annihilation, 
creation and conservation processes: 

A k (t)e(f)= f f k ( S )dse(f), 
Jo 

(e(g)\Al(t)e(f )) = f ^(s) ds (e(g)\e(f)), 
Jo 

(e(g)\A kl (t)e(f)) = f M /, (s) ds (e(g)\e(f)). 
Jo 

The annihilation and creation processes satisfy the canonical commutation rules 
(CCR); formally, [A k (t) , A] (s)] = t A s, [A k (t), A^s)} = 0, [A{(t) , A\ (s)] - 0. 

Let H, R k , S k i, k, I = 1, . . . , d, be bounded operators on Tt such that H* = H 
and S* k S jt - J2 3 SkjStj - 5 k i. We set also K := -iH - ± £ fe R* k R k . Then, 
the quantum stochastic differential equation 17 

d(7 W = IE Rk dA l^ + E ( Skl - S ki)dA kl (t) 

^ k kl 

-J2R* k SkidMt)+Kdt\u(t), (1) 

kl ' 

with the initial condition U(0) = 1, has a unique solution, which is a strongly 
continuous family of unitary operators on H <£> T, representing the system-field 
dynamics in the interaction picture with respect to the free field evolution. 
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1.2 The reduced dynamics of the system 

The states of a quantum system are represented by statistical operators, positive 
trace-class operators with trace one; let us denote by S(7i) the set of statistical 
operators on H. As initial state of the composed system "system S plus fields" 
we take p <8> f?r(/), where p G S{TL) is generic and Qr{f) is a coherent state, 
Qr(f) '■= |e(/))(e(/)|. One of the main properties of the Hudson-Parthasarathy 
equation is that, with such an initial state, the reduced dynamics of system £ 
obeys a quantum master equation. 16, 17 Indeed, we get 

^Th = £{t)fa], Vf.= Tr r {U(t)(p® Qr (f))U(ty}, (2) 
where the Liouville operator C(t) turns out to be given by 

m\p] = [K-YsKSkiMt)) p+pl^-j^ms^Rk 

\ kl / \ kj j 

+ E - E p (5 - s k*w) ii/wii V (3) 

A particularly important case is Ski — Ski, when £(t) reduces to 

k k 

+ UkPRl - Wk*kP - \pR* k Rk) ■ (4) 

It is useful to introduce also the evolution operator from s to t by 

^T(t,i) = £(i)oT(t,«), T(«,a) = l. (5) 
With this notation we have rj t = T(t,0)[p]. 

1.3 The field observables 

The key point of the theory of continual measurements is to consider field ob- 
servables represented by time dependent, commuting families of sclfadjoint op- 
erators in the Heisenberg picture. 16 Being commuting at different times, these 
observables represent outputs produced at different times which can be obtained 
in the same experiment. Here we present a very special case of family of ob- 
servables, a field quadrature. Let us start by introducing the operators 



Q(t- 0,i/)= f e-^ s+ ^dA\(s)+ f e i ("" + *><L4i(s), t > 0; 
Jo Jo 



(6) 
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i3 6 (— 7r, 7r] and v > are fixed. The operators Q(i; v) are selfadjoint (they 
are essentially selfadjoint on the linear span of the exponential vectors). By 
using CCR's, one can check that they commute: [Q(t; v), Q(s; d, v)] = 
(better: the unitary groups generated by Q{t\ , d,v) and Q(s;i?, v) commute). 
The operators |B]| have to be interpreted as linear combinations of the formal 
increments cL4.J(s), &A\(s) which represent field operators evolving with the 
free-field dynamics; therefore, they have to be intended as operators in the 
interaction picture. The important point is that these operators commute for 
different times also in the Heisenberg picture, because 

Q™\t-d,v) := U{t)*Q{t;&,v)U{t) = U(T)*Q(t;d,v)U(T), VT > t; (7) 

this is due to the factorization properties of the Fock space and to the prop- 
erties of the solution of the Hudson-Parthasarathy equation. These "output" 
quadratures are our observables. They regard those bosons in "field 1" which 
eventually have interacted with S between time and time t. Commuting 
selfadjoint operators can be jointly diagonalized and the usual postulates of 
quantum mechanics give the probabilities for the joint measurement of the ob- 
servables represented by the selfadjoint operators Q oat (t; i?, v), t > 0. Let us 
stress that operators of type © with different angles and frequencies represent 
incompatible observables, because they do not commute but satisfy 

rnr , a \ r\t ± \\ 4isin(^- M ))sin(g-0+^Q,- M )) 
[Q(t; 9, v),Q{s; 0, fi)] = ^ '- ^ ^ '- . 

v — fl 

When "field 1" represents the electromagnetic field, a physical realization 
of a measurement of the observables ([7]) is implemented by what is called an 
heterodyne/homodyne scheme. The light emitted by the system in the "channel 
1" interferes with an intense coherent monochromatic laser beam of frequency 
v. The mathematical description of the apparatus is given in Section 3.5 of Rcf. 
16. 



1.4 Characteristic operator, probabilities and moments 

The commuting selfadjoint operators ^ have a joint pvm (projection valued 
measure) which gives the distribution of probability for the measurement. Any- 
way, at least the finite-dimensional distributions of the output can be obtained 
via an explicit and easier object, the characteristic operator <3? t (fc;i9, v), a kind 
of Fourier transform of this pvm. For any real test function k 6 i°°(R + ) and 
any time t>0we define the unitary Weyl operator 

$ t (fe;i?,i/) = cx p|i^ fc(s)dQ(s;tf,2/)j 

= ex p{i^ k{s)e-'^ s+ ^AA\{s) fc(s)e i ( l/s+,? )dA 1 ( S )|. (8) 
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Then, there exists a measurable space (O, J 7 ), a pvm £^ (acting on T) with value 
space (fi,^ 7 ), a family of real valued measurable functions {X(t; •) , t > 0} on 
fi, such that X(0;w) = 0, and, for any choice of n, = to < ti < • • • < t n < t, 



$ t (fc;i?,i/) = cxpl i^2nj[Q(tj;tf,v) - <3(t,_i; i/) 

= / fi ^{'^ Kj W) " W) ] }^ (dw) ' (9) 

where k(s) = Y^j=i l(t 3 -i,t,)( s ) K j- Let us stress that the pvm depends on 
the observables and, so, on the parameters z? and v, while the choices of the 
trajectory space (the measurable space (Q,.F)) and of the process X(t) are 
independent of ■& and v. 

Then, we introduce the characteristic functional 



t 

,out / 



$*(*;; 0,1/) =Tr|cxp|i j fc(s) dQ out (s; 0, i/)|p ® g r (f) 

= Tr {$ t (fc; 0, i/)I7(t) (p ® ??(/)) <7(t)*} . (10) 

All the finite-dimensional probabilities of the increments of the process X(t) are 
determined by 

$t(fc;0,i/) = ( ll'-" " 

x P^[AX(i ,*i) G dxi,...,AA'(t n _i,t n ) e <Lc„] , (11) 

where we have introduced the test function k(s) — Y^j=i l(t,_i,t,)( s ) K j) with 
= t < *i < • • • < t n < t, Kj e R. 

The fact that the theory gives in a simple direct way the distribution for the 
increments of the process X(t), rather than its finite-dimensional distributions, 
is related also to the interpretation: the output X(t) actually is obtained by a 
continuous observation of the generalized process I(t) — dX(t)/dt followed by 
post-measurement processing. 

Starting from the characteristic functional it is possible to obtain the mo- 
ments of the output process I(t) and to express them by means of quantities 
concerning only system S. 16 Let us denote by ¥f p ,v the expectation with respect 
to Pp'"; for the first two moments we obtain the expressions 

E^[I(t)}=Tr{(Z(t) + Z(ty) Vt }, (12a) 



E^[I(t)I(s)}=S(t-s) 

+ Tr {(Z(t 2 ) + Z(t 2 )*) T(t 2 , ti) [Z(t x ) mi + Vt.Zih)*}} , (12b) 
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where t\ = t A s, t<i = t V s and 

Z(t):= e^ t+ ") + £ Si*/* (i) j . (12c) 

2 The spectrum of the output 

2.1 The spectrum of a stationary process 

In the classical theory of stochastic processes, the spectrum is related to the 
Fourier transform of the autocorrelation function. Let Y be a stationary real 
stochastic process with finite moments; then, the mean is independent of time 
E[Y"(i)] = E[y(0)] =: my, Vt € R, and the second moment is invariant under 
time translations 

E[Y(t)Y(s)]=E[Y(t- s)Y(0)]=: R Y (t- s), Vi,seR. (13) 

The function Ry(t), t S R, is called the autocorrelation function of the process. 
Obviously, we have Cov [Y(t),Y(s)] = Ry(t— s) - my. 

The spectrum of the stationary stochastic process Y is the Fourier transform 
of its autocorrelation function: 

/+oo 
e iMr i?y(r)dT. (14) 

This formula has to be intended in the sense of distributions. For instance, if 
Cov [Y(t), y(0)] G L 1 ^), we can write 

r + co 

S Y {n) := 2Trm^5(fi) + / e^ r Cov [F(r), Y (0)] dr . (15) 



By the properties of the covariance, the function Cov [Y (r), 1^(0)] is positive 
definite and, by the properties of positive definite functions, this implies 

e if * T Cov [r(r),y(0)]dr > 0; 

J — oo 

then, also Sy{[i) > 0. 

By using the stationarity and some tricks on multiple integrals, one can 
check that an alternative expression of the spectrum is 



S y (jm)= Hm i E 

T^+oo 1 



T 

e^Yit) dt 



(16) 



The advantage now is that positivity appears explicitly and only positive times 
are involved. Expression (|16[) can be generalized also to processes which are 
stationary only in some asymptotic sense and to singular processes as our I(t). 



6 



2.2 The spectrum of the output in a finite time horizon 

Let us consider our output I{t) = dX(t)/dt under the physical probability ¥® p M . 
We call "spectrum up to time T" of I{t) the quantity 



T 



e 1 "* dX(t) 



(17) 



When the limit T — > +00 exists, we can speak of spectrum of the output, but 
this existence depends on the specific properties of the concrete model. 

By writing the second moment defining the spectrum as the square of the 
mean plus the variance, the spectrum splits in an elastic or coherent part and 
in an inelastic or incoherent one: 



Mm; v) = S${p\ 0, v) + S™ 1 ^, d, v), 











[ e 1 ^* dX(t) 






Jo 





(18) 
(19) 



5r 1 (/i;0,z,) = ivar^ 



cos jit dX(t) 



rp p 



sin jit dX(t) 



, (20) 



Let us note that 

SfJ(/x; 0, u) = 5|}(-/x; 0, u), S™ l (ii; 0, 1/) = ( 21 ) 



By using the expressions (| 1 2[) for the first two moments we get the spectrum 
in a form which involves only system operators: 



&*Tr {(Z(t) + Z(t)*)r) t }dt 



(22a) 



S^ c \ii] 0, 1/) = 1 + - / dt / ds cos fj,(t - s) 
T Jo Jo 



x Tr { (z(t) + T(t, s) + ffrZ(«)*] } , (22b) 

Z(t) = Z(t) - Tr {£(*)??*} . (22c) 
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2.3 Properties of the spectrum and the Heisenberg uncer- 
tainty relations 

Equations (|22[) give the spectrum in terms of the reduced description of system 
S (the fields are traced out); this is useful for concrete computations. But the 
general properties of the spectrum are more easily obtained by working with 
the fields; so, here we trace out first system S. Let us define the reduced field 
state 

n T (/) := Ttn {U(T)(p® Qr (f))U(T)*} (23) 
and the field operators 



V-t Jo 

Qr(« 0, v) = Qt(m; i?, v) - Tr {U T (f)Q T (^; 0, v)} . 



(24a) 
(24b) 



Let us stress that Qr(/x; i 7 ) commutes with its adjoint and that Qt(m; ^)* = 
Qt{~^; J?, By using Eqs. (jTO]) and (fTTj) and taking first the trace over Ti., we 
get 

St&I 0, i/) = Tr {n T (/)Q T ( M ; 7?, ^)*Q t (m; ^, *)} > 0, (25a) 
i?, i/) = |Tr {n T (/)Q T ( M ; 0, ^)}| 2 > 0, (25b) 

S^ nel (/i; ??, i/) = Tr {n T (f)Q T (n; 0, ^)*Qt(m; i?, f)} > 0. (25c) 

To elaborate the previous expressions it is useful to introduce annihilation 
and creation operators for bosonic modes, which are only approximately orthog- 
onal for finite T: 



e i " t dA 1 (t + T/2), 



(26a) 



[a T (u),a T (u')] = [a^(w),a^(w')] = 
[a T (w),a^(w')] : 



1 for lu' = u), 



for u' 7^ u). 



Then, we have easily 



(26b) 
(26c) 

(27) 



5t(m; 1?, = 1 + Tr|n T (/)^a^(^ + ^)a T {v + fj) + a^ T (i> — fi)a T (v - (Jt) 

+ e~ m a! r {v + n)a\,(v - /i) + e 2W a T (^ - ^)a T {v + /*))}> ( 28a ) 

0,") = le^TrlnT^aT^ + ^j + e-^TrlnT^^-/!)}! 2 . (28b) 
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Theorem 1. Independently of the system state p, of the field state gr{f) and 
of the Hudson- Parthasarathy evolution U , for every {> and v we have the two 
bounds 

\ (S™\^ 0, v) + S™Hp; ± §, *)) > 1, (29) 

s^ nc V; 0, v)S't c \h; ± f , v) > l. (30) 

Proof. The first bound comes easily from 

S™\v 0, i/) = 5 t (m; i?, «/) - S^(w 0, f) 



and Eqs. H2HJI. 

To prove the second bound, let us introduce the operator 

B T (u>) :=a T (i/-w)-Tr{n T (/)a T (j/-w)}, (31a) 
which satisfy the CCR 

[Bt(u),bUu)] = 1, [Bt{u), B t {lo)] = {bU"), 4M] = 0- (31b) 
Then, we can write 

S*? el (/x; 0, „) - Tr{ (e^4(-^) + e^Brfc)) O r (/) 

x ( e - M 4M + e M B T (-p))}. 

The usual tricks to derive the Heisenberg-Scrodinger-Robertson uncertainty 
relations can be generalized also to non-selfadjoint operators. 3, 18 For any choice 
of the state g and of the operators X\, X 2 (with finite second moments with 
respect to g) the 2x2 matrix with elements Tr {Xi£>X?} is positive definite 
and, in particular, its determinant is not negative. Then, we have 

Tr{X 1 gX*}Tr{X 2 gX*} > \Tr{X lQ X*}\ 2 

> \ImTT{X ig X;}\ 2 = \ \St{ Q {X\X x - X*X 2 )}\ 2 . 

By taking g = IL T (f), X x = B^(p) + e i€ B T (-fi), 
X 2 = T i (e- w Bl(ji) - e^B T {-fj^ , we get 



SF A (j*;0,u)SF*( M ;0±%,u) 
> 



1 + Tr {n T (/) (4( M )B T ( M ) - 4(-/i)S T (-/x)) } 
But we can change \i in — \i and we have also 

s™\»; 0, v)s™\^ ± f , = sr\-»; 0, ^)4 nel (-M; ± f , 

1 + Tr |n T (/) (4(- M ) J B T (-/i) - 4( M )£ t (m)) } 



> 
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The two inequalities together give 



2 



> 1 



Tr{n r (/)(4(/x) J B T ( At )-4(-/i) J B T (- At ))}|) >1, (32) 

which is what we wanted. □ 

Ref. 18 introduces a class of operators for the electromagnetic field, called 
two-mode quadrature-phase amplitudes, which have the structure (|27[) of our 
operators Qt{h\ v). Anyway only two modes are involved, as if we fixed /i and 
v. Let us denote here those operators by Q pa . The paper explicitly constructs 
a class of quasi-free (or Gaussian) field states g sq for which Tr {g S qQp a Q P a} — 
|Tr {&qQpa}| 2 < 1. Such states are called two-mode squeezed states. More 
generally, one speaks of squeezed light if, at least in a region of the /i line, for 
some $ one has S^ cl (fi; v) < 1. If this happens, the Heisenberg-type relation 
(|3"0|) says that necessarily S™ cl (n; i? + v) > 1 in such a way that the product 
is bigger than one. 



3 Squeezing of the fluorescence light of a two- 
level atom 

Let us take as system S a two- level atom, which means TC = C 2 , H — ^ a z ; 
loq > is the resonance frequency of the atom. We denote by cr_ and <r+ 
the lowering and rising operators and by o~ x = o~_ + cr + , a y = i(<r_ — 0+), 
ct z = cr + cr_ — ct_(t + the Pauli matrices; we set also a® = e 1 <r_ + e - "' cr + . We 
stimulate the atom with a coherent monochromatic laser and consider homodyne 
detection of the fluorescence light. The quantum fields T model the whole 
environment. The electromagnetic field is split in two fields, according to the 
direction of propagation: one field for the photons in the forward direction 
(k = 2), that of the stimulating laser and of the lost light, one field for the 
photons collected to the detector (k = 1). Assume that the interaction with 
the atom is dominated by absorption/emission and that the direct scattering is 
negligible: 

Ski — Ski , i?i = ^/tp <t_ , R 2 = V7(l -p)a- . 

The coefficient 7 > is the natural line-width of the atom, p is the fraction of 
fluorescence light which reaches the detector and 1 — p is the fraction of lost 
light (0 < p < l). 10 ' 15, 16 ' 19 We introduce also the interaction with a thermal 
bath, 

i?3 = ^Jyho'- , i?4 = \J "fna + , n > 0, 
and a term responsible of dephasing (or decoherence), 

R5 = V ikd o~ z , k d >0. 
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To represent a coherent monochromatic laser of frequency to > 0, we take 



fk(t) 



Jk2 



f l[o,T]W; T is a time larger than any other time in 
-hoo is taken in all the physical quantities. The 



the theory and the limit T 
quantity O > is called Rabi frequency and Aw = loq — lu is called detuning. 
The squeezing in the fluorescence light is revealed by homodyne detection, which 
needs to maintain phase coherence between the laser stimulating the atom and 
the laser in the detection apparatus which determines the observables Q(t; v); 
this in particular means that necessarily we must take v = lu. 

The limit T — > +00 can be taken in Eqs. ([22]) and it is independent of the 
atomic initial state. 19 The result is 



S c \fi;tf):= lim 



(33) 



where 



S* mo V; 1?) := lim S^ no V; 0, w) = 1 + 2 7 p 

T — ^+oo 



.4 



A 2 + fi 2 



(34) 



t = Tr 



(i ?9 — i i9 

e er_p cq + p cq e c 



Peq 



A = 



2k A 









1 cost? 


Pcq] Pcq) <? 


1 




sin$ 








V 










x cq = -jA~ 


1 


(!)■ 












Aw 




° \ 




+ n + 2k d ) 













7(1 + 2n), 



Examples of inelastic spectra are plotted for 7 = 1, n = kd. = 0, p = 4/5, 
and two different values of Aw. The Rabi frequency Q and 9 are chosen in 
both cases to get good visible minima of S inel below 1. Thus in this case the 





Figure 1: S* inol (/i; i9) with Aw 
3.5, = 3.7021. 



Figure 2: S inel (iJ,; 1?) with Aw = 0, 
= 0.2976. 
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analysis of the homodyne spectrum reveals the squeezing of the detected light. 
Also complementary spectra are shown to verify Theorem 2.1. One could also 
compare the homodyne spectrum with and without n and fc<j, thus verifying 
that the squeezing is very sensitive to any small perturbation. 
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